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A new method is introduced to study three-body clusters. Triangular configurations with T>3h 
point-group symmetry are analyzed. The spectrum, transition form factors and B(EX) values of 12 C 
are investigated. It is concluded that the low-lying spectrum of 12 C can be described by three alpha- 
particles at the vertices of an equilateral triangle, but not as a rigid structure. Large rotation-vibration 
interactions, Coriolis forces and vibration-vibration interactions are needed. Other interpretations, 
such as the harmonic oscillator and a soft deformed oscillator with SO(6) hyperspherical symmetry, 
appear to be excluded by electron scattering data. 
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Abstract 



1 Introduction 



X 

The Cluster Model of nuclei was introduced at the very beginning of Nuclear Physics as a natural 
consequence of the molecular viewpoint in which aggregates are formed by combining constituent units 
into composite objects. In particular, the alpha-particle model represented, for many years, a viable 
model of light nuclei |Q. With the successes of the Nuclear Shell Model, the Cluster Model lost most of 
its appeal and the interest shifted in 'deriving' the alpha-particle model from the Nuclear Shell Model 
[pj.A large literature exists in this direction, especially from the Japanese school 1|, 0, ||. In recent years, 
clustering in nuclei has found renewed interest, since it appears that many new isotopes of light nuclei 
far from stability possess a cluster structure (not necessarely alpha). In view of this renewed interest, we 
have readdressed the question of clustering in nuclei with two purposes in mind. First, to analyze cluster 
configurations at the 'classical' level in somewhat more detail than previously done, deriving explicit 
expressions for the density distribution and its multipole moments, and second, and most importantly, to 
introduce a formalism that allows one to study in a relatively simple fashion the actual situation which 
occurs in nuclei, where clustering is not a rigid-like molecular structure, but a rather soft (liquid-like) 
structure. The new approach is based on the bosonic quantization of many-body problems introduced 
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by one of us Q. In this method, a quantum mechanical problem with v space degrees of freedom is 
quantized in terms of elements of the Lie algebra U{y + 1). This method has been used in a variety of 
problems, ranging from Nuclear Physics to Molecular Physics, from Hadronic Physics to Polymer Physics 
and has proven to be very useful in analyzing experimental data. In a series of articles beginning with 
the present one, we intend to exploit the algebraic method for Cluster Physics. Two-body clusters, where 
the algebra is C/(4), have been discussed previously Q|. Here we begin with three-body clusters, where 
the algebra is Z7(7), and apply the method explicitly to the study of three-alpha configurations in nuclei, 
in particular to 12 C. Many authors have suggested that this nucleus in its ground state is composed 
of three alpha-particles at the vertices of an equilateral triangle. In this article we study, both at the 
classical and at the quantum level, a configuration of three-particles at the vertices of an equilateral 
triangle (point-group symmetry T> 3 h) and analyze the extent to which the experimental data support 
this configuration. A preliminary account of this part of our work has appeared g, 0]. In addition to 
the study of properties of triangular configurations we also discuss other possible arrangements of three 
particles and their motion, in particular those corresponding to a six-dimensional spherical oscillator and 
to a six-dimensional deformed oscillator. Altough these configurations do not appear to describe well the 
experimental data in 12 C the formulas we derive may be of interest in other situations involving three 
particles. 

2 Classical treatment 

The treatment of a triangular configuration in classical mechanics is very well known from molecular 
physics. However, this treatment is usualy done in terms of point-like constituents. Although the gener- 
alization to composite constituents is straighforward, in this introductory section we present some new 
and more detailed results. 

2.1 Energies 

We consider three identical particles of mass Am/ 3 at the vertices of an equilateral triangle (point-group 
symmetry T>^h) (see Fig. [l]). The origin is placed at the center of mass and the three particles span the xz 
plane. The distance from the center is /3. The spherical coordinates of the three particles r"i = (r^, 9i,4>i) 
are then (/3,0, 0), (/3, 2tt/3, it) and (/3, 27r/3,0) for i = 1,2 and 3, respectively. For a mass distribution 
p(r), rotational energies can be computed by evaluating the moments of inertia 

T 
T 

y 

The rotational energy is simply 

2.1.1 Point-like masses 

For a point-like mass distribution 

p(f) 

i=l 



(y 2 + z 2 )p(f) dr 
(z 2 + x 2 )p(r) dr 
(x 2 + y 2 )p{r) dr 



(1) 



1 t2 , 1 t2 , 1 r2 

2x: L * + 2z:^ + 2x: L2 
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a simple calculation gives 



= -Am.0 2 , 
2 

= Amp 2 , 



and thus 



Am(3 2 



L 2 



(4) 
(5) 



where L is the angular momentum and K its projection on the three-fold y-axis. 
2.1.2 Extended mass distribution 

Consider the case in which instead of point masses there are, at the vertices of the triangle, extended 
distributions. Since we have in mind applications to the alpha-particle model, we consider here the case 
in which the distribution is Gaussian 



p(r) 



Am (a\3/2^ r 5 



i=i 



(G) 



It is straightforward to compute the moments of inertia 



I. 



\Am/3 2 ( I 
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The rotational energy is now given by 

Erot = — 



Am/3 2 ( 1 + - 



L 2 — 



af3 2 



2(1 + a/3 2 



-K 1 



(7) 



(8) 



For typical values of a and (3, as for example those extracted from the form factors of 12 C (discussed in 
the following sections), a — 0.52 fm -2 and (3 = 1.74 fm, we find that the coefficient in front of K 2 reduces 
from —0.50 to —0.30. We note in general that the rotational energy of a symmetric top is given by 



E r , 



= AL 2 + BK 2 



(9) 



The top is called oblate if B < 0, prolate if B > and spherical if B = 0. Three particles on a plane are 
always oblate. The extended distribution tends to make the top less oblate. 

Vibrational energies can be computed by tying the particles with strings and evaluating the eigenfre- 
quencies . Denoting by q\ , q2 and q$ the displacements of the sides of the triangle and assuming a 
potential energy of the type 

\C (qj + q 2 +q 2 )+D { qi q 2 + q 2 q 3 + q m ) , (10) 



V 



one can obtain the frequencies of the normal modes as 



to i 



^2 



9{C + 2D) 



Am 



9(C-D) 
2Am 



(11) 



The vibration u)% is singly degenerate, while the vibration 0J2 is doubly degenerate. 
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2.2 Transition probabilities 

Consider now the case in which the three particles have a charge Ze/3 (the case in which the particles 
have a magnetic moment will not be discussed here since we have in mind the alpha-particle model, 
where the constituents have no magnetic moment). The total charge distribution can be expanded into 
multipoles 

P (r) = Y, A im{r)Y lm {e,cj>) , (12) 

with 



A lm {r) = J Y? m (d,<t>)p(f)dcoB9<fy . (13) 
Instead of the charge distribution, one can consider the form factor, given by its Fourier transform 

F(q) = ±J e^p(r)df = J2 F ^n(q)Y em (q) , (14) 



trn 



Fe m (q) = -7T I r 2 A hn (r)i l je(qr)dr . (15) 



with 

(n\ = 

Ze 

The multipole moments of the charge configuration are 

Qe m = Jr l Y e * m (9, 4>) p(f) dr = / A £m (r) r l + 2 dr 

= £f(2l+l)!!lim%M. (16) 
47T <j-»o % t q t 

From these multipole moments, the transition probabilities can be calculated classically |fll|| . The tran- 
sition probability per unit time is 

T(E£) = 8nc^- — 5 k 2£+1 B(E£) , (17) 



where 



he e[{2£ + iy.\] 



B(El) = J2 QlmQtrn ■ (18) 



2.2.1 Point charges 

For a point-like charge distribution 



the form factor is given by 



P(f) = ?E^-^), (19) 



3 i 



F(q) = 5 (20) 



i=l 
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The corresponding B(E£) values are 



which gives 



B(E£) = — ) /} 



Ze\ 2 n 2£+l 



47T 



3 + 6P/(--) 



(21) 



B(E2) = (Zefl\0\ 
B(E3) = (ZefHf3\ 

The dipole radiation B(E1) vanishes because of the spatial symmetry of the charge distribution. 



(22) 



2.2.2 Extended charge distributions 

Next we consider the case in which instead of point charges there are, at the vertices of the triangle, 
extended distributions 



exp 



-a (r — fi)' 



Ze /a\ 3 / 2 
~3 



3 /2 

(^) e- Q(r2+ ' 32) 47r^^ iA (2 a /3r)r A (0,^.y;(^,^) , 



(23) 



i=l A=0 



where u(a;) = jx{ix)/i x is the modified spherical Bessel function. The corresponding form factor is that 
of the point-like charge distribution multiplied by an exponential 



F(q) 



(24) 



i=i 



The B(Et) values of the extended charge distribution are the same as those of the point-like configuration. 



3 Quantum treatment 



The quantum treatment of an identical three-body cluster can be done in several ways. The most 
straightforward is to introduce relative Jacobi coordinates 



P = (n-r 2 )/V2, 

A = (r 1 +f 2 -2f 3 )/V6 



(25) 



write down a Hamiltonian in terms of these coordinates and their canonically conjugate momenta p p and 
p\, and solve the Schrodingcr equation 



1 



V(p,A) 



Ei,{p,\) 



(26) 



to obtain the energy eigenvalues and eigenvectors. 
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3.1 Energies 

It is convenient to make a change of variables to coordinates (r, £, f2 p , Q\) which are related to the 
Jacobi coordinates (p, Q, p , A, by 



p = rsin£ , A = rcos£ 

In the new coordinates, the kinetic energy is given by fl^, 0| 



1 



d 



2mr 5 dr 



,9_ 
dr 



1 



2mr 2 



—2 -4cot2£— + : 2 
9^ 3£ sin 2 £ 



cos 2 £ 



(27) 



(28) 



In this section, we consider three different potentials, namely the harmonic oscillator, the deformed 
oscillator and the oblate symmetric top. For potentials that only depend on the hyperradius r, the 
Schrodinger equation can be solved by separation of variables into an angular and a radial equation 



_ 4co t2£— + -^Vr + ^Vr 

dt; 2 <9£ sin ^ cos 2 ^ 



1 1 d 

2m r 5 dr 



■A 

dr 



o-(er + 4) 



2mr 2 



For the six-dimensional harmonic oscillator 



V(r) 



V(r) 



:Cr 2 



tp(r) 



the energy spectrum can be obtained exactly 

E(n) = e {n + 3) , 



Eip(r) . (29) 

(30) 
(31) 



with n — 0, 1, . . . , and e = y/ C/m. The allowed values of a are a = n, n — 2, . . . , 1 or for n odd or 
even, respectively. 

For the six-dimensional displaced (or deformed) oscillator 



V(r) 



;C(r-r ) 2 , 



(32) 



the energy eigenvalues cannot be obtained exactly. An approximate solution, valid in the limit of small 
oscillations around the equilibrium value rg is 



E(v,a) 



e {v 



1 



2mr 2 



cr(er + 4) + 



15 



(33) 



with e = a/C/to. The first term gives a harmonic vibrational spectrum with v = 0, 1, . . . , and the second 
term gives the rotational spectrum with a = 0, 1, . . . . 

In general, the potential is not invariant under six-dimensional rotations as in the previous two 
examples, but only under rotations in three dimensions. One can introduce the total angular momentum 
I = L p + L x , and make a transformation of variables to the relative angle 26 between the two vectors 



cos 26 = p ■ A , 

and the Euler angles Q of the intrinsic frame defined by |l4| 

p x A - p — A 



1 



sin 26 



2sin( 



p + X 
2cos6» 



(34) 



(35) 
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In the new variables, the angular momentum dependent terms in the kinetic energy are given by 



L 2 p (Q p ) t L\(n x ) 
cos 2 £ 



1 



sin 2 £ 



4 sin £ cos 2 £ 



d 2 d 
tan 2 9 if + cot 2 if - — r - 2 cot 26>— 



cos2£ 



cot 6 Inl-i — tan( 



2 sin £ cos 2 £ 

The potential only depends on the intrinsic variables r, £ and 0. 
the potential has sharp minima in r, £ and 



(36) 



An interesting situation occurs when 



\c{r 



9o) 2 



(37) 



In the limit of small oscillations around ro, £o — 7r /4 and 9q = tt/4, rotations and vibrations decouple, 
and the set of resulting differential equations can be solved in closed form 



E(vl,V2a,V2b,I,K) = 61 (Ul + i) + e 2a («2a + ^) + £2fc («2fc + Jj) 



1 



-if 2 



(38) 



with ei = yJCjm, e 2a = \J A/ mr* and e 2 {, = \JBjmr^. Here if is the projection of the angular 
momentum I on the symmetry axis 1. 

For three identical masses, the potential has to be invariant under their permutation, i.e. the coeffi- 
cients A and B are equal. As a consequence, the vibrations i> 2a and v%\, form a doubly degenerate vibration, 
and it is convenient to use instead of V2a and V2b the quantum numbers v 2 2 where £2 — v 2 , V2 — 2, . . . , 1 
or for V2 — « 2 a + V2b odd or even. The vibrational energies can then be written as 

E vih (v 1 ,v 2 ) = ei (wi + -) + e 2 (v 2 + 1) . (39) 



The rotational energies 



E rot (I,K) 



HI 



(40) 



are the same as derived classically for a triangular distribution of three identical point-like masses of 
Eq. (H) with I 2 replaced by 1(1 + 1). The spectrum is that of an oblate symmetric top. 

The triangular configuration with three identical masses has discrete symmetry T>3h.. Among the 
quantum numbers, it is convenient to add the parity P and the transformation properties of the states 
under T>3h- Since T>3h ~ x P, the transformation properties under T>3h are labeled by parity and 
the representations of V3 ~ S3. V3 has two one-dimensional representations, usually denoted by A\ 
and A 2 , and a two-dimensional representation denoted by E. The corresponding Young tableau of the 
permutation group S3 are 



t = A 1 



(41) 



t = E 



(42) 



t = An 



(43) 



For three identical particles, the wave functions must transform as the symmetric representations A\ of 

lh ~ s,. 
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3.2 Transition probabilities 

Transition probabilities, charge radii and other electromagnetic properties of interest can be obtained 
from the transition form factors. For electric transitions (the only ones discussed here) the form factors 
are the matrix elements of the Fourier transform of the charge distribution 

F(i -> /; = J dve 1 ^ ( 7/ , L f , M f | p(f) | 7i , Li, M z ) , (44) 

where 7 denotes the quantum numbers needed to classify the states uniquely. For the point-like charge 
configuration of Eq. (|l9|), the transition form factor reduces to 

F(i —> f;q) = ZeY,V%;U$ni^fiQ)V%& i {-Q), (45) 

M 

with 

F(i^f;Q) = (T/.i/.MIe-W^-^.Li.M) . (46) 

Here we have used the permutation symmetry of the initial and final wave functions, a transformation to 
Jacobi coordinates, and an integration over the center-of-mass coordinate. In Table [l] we show some form 
factors for the transition from the ground state to an excited state !F(0^ — ► Lf ; q) . The form factors 
for harmonic oscillator show an exponential dependence on the momentum transfer jl^, , whereas for 
the deformed oscillator and the oblate top they behave as cylindrical |l7| and spherical Bessel functions 
H [H], HU , respectively. The form factors associated with the extended charge distribution of Eq. (|23| ) can 
be obtained, as in the classical case, by multiplying the results from Table [j]by a Gaussian exp(— q 2 /4a), 
which represents the form factor of the a particle. 

The transition probability B(EL) can be extracted from the form factors in the long wavelength limit 



B{EL:i^f) = [(2Z/ + 1)!!] ' lim V \F(i -> /; q)\ 2 /q 
1 ' J ' 442^ + 1) r+oJ-^J V n 1 



2L 



w^sffinXw-r.tf/r. (47) 



In the case of the oblate symmetric top, we find for the rotational excitations of the ground state the 
same results as for the multipole radiation in the classical case (see Eq. (p2|)) 



B(E2;0+ - 


.2+) 




B{E3;0 + - 






B(E4;0 + - 


-4+) 





(48) 

The ground state rotational band of the oblate top does not have a symmetric state with L p — 1~, which 
is in agreement with the absence of dipole radiation in the classical treatment. 

The ground state charge distribution can be extracted from the elastic form factor by taking the 
Fourier transform 

p(r) = — L, f dqF(0+ ^0+;q)e-^ 



(2tt) 3 
Ze 
2^2 



q 2 dqj (qr)F(Oi ^0+;q) ■ (49) 
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For the case of the oblate top, the ground state charge distribution can be derived as 

3/2 i 



p(r) 



Ze 



(!)' 



-a(p-r) 2 



-a(0+ry 



4a (3r 

The charge radius can be obtained from the slope of the elastic form factor in the origin 

1/2 

(r 2 ) 1/2 



-6 



dF(0+->0+;<z) 



dq 2 



q=0 



(50) 



(51) 



In Table [l] we have introduced a scale parameter (3 in the form factors, such that the charge radius for 
the harmonic oscillator, the deformed oscillator and the oblate top is given by 



,2U/2 = 



2a 



d 2 



(52) 



4 Algebraic treatment 

In this section, we discuss an algebraic treatment of three-cluster systems, in which the eigenvalue problem 
is solved by matrix diagonalization instead of by solving a set of differential equations. In the algebraic 
cluster model (ACM), the method of bosonic quantization is used. This method consists in quantizing 
the Jacobi coordinates and momenta with boson operators and adding an additional scalar boson Jig] 



u p,m J u x, m 1 A 



t 



: m = 0,±1) 



(a = l,...,7) 



(53) 



The role of the additional scalar boson is to compactify the space thus making calculations of matrix 
elements simpler. The set of 49 bilinear products of creation and annihilation operators spans the Lie 
algebra of U(7) 



clc/3 



(a,/3=l,...,7) 



(54) 



All operators are expanded into elements of this algebra. 

The most generic Hamiltonian which includes up to two-body terms and is invariant under the permu- 
tation group S3, contains 10 parameters that describe the relative motion of a system of three identical 
clusters. It can be written as [115] 



H 



e s+s - e 1 (fit ■ b p + b\ ■ b x ) + u (sVss) - u x s^&J • b p + b[ ■ b x )s 
+v [(bl -b\ + b{- b\)s~s + »V(g p . ~ bp + ~ bx . b x ) 

+ J2 Wl ( b l xb l + b l x b l) (l) ■ % x ~ h P + ^ x 6 A ) (0 



1=0,2 



+ 



Q [(bl x fit _ 6 t x fit)W . (b p xb p -b x x 6 A )« + 4 (^ x &t)« . (i x x ~ bp) 



1=0,2 



+cx (bl x b{)^ ■ (b x x b p )W , 



(55) 



with b n 



( — 1 ) 1 m bp,- m , b\ jn = (— l) 1 m b X - m and s = s. Here the dots indicate scalar products 



and the crosses tensor products with respect to the rotation group. By construction, the Hamiltonian 
commutes with the number operator 



N 



sU-bi-L-bX-bx 



(56) 
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and therefore conserves the total number of bosons N — n s + n p + n\. In addition to N, the eigenfunctions 
have good angular momentum L, parity P, and permutation symmetry t. The three symmetry classes 
of the S3 permutation group are characterized by the irreducible representations: t = S for the one- 
dimensional symmetric representation, t = A for the one-dimensional antisymmetric representation, and 
t = M for the two-dimensional mixed symmetry representation. 



4.1 Dynamic symmetries 

The 5*3 invariant U (7) mass operator of Eq. (|55|) has a rich algebraic structure. It is of general interest 
to study limiting situations, in which the energy spectra can be obtained in closed form. These special 
cases correspond to dynamic symmetries of the Hamiltonian. Dynamic symmetries arise, when (i) the 
Hamiltonian has an algebraic structure G, and (ii) it can be expressed in terms of Casimir invariants of a 
chain of subalgebras of G only. Here we consider two different dynamic symmetries, corresponding to the 
six-dimensional spherical and deformed oscillator cases discussed in the Schrodinger picture in Sect. 3. 

For vq = in Eq. (|5^), we recover the six-dimensional oscillator model, since there is no coupling 
between different harmonic oscillator shells. The oscillator is harmonic if all terms except eo and ei are 
set to zero, otherwise it is anharmonic. This dynamic symmetry corresponds to the reduction 

U(7)DU(6),D--- (57) 

The dots here indicate further algebras needed to classify the states uniquely. Since our purpose here 
is to present a global analysis of the situation, we do not dwell on the group theoretical aspects of this 
problem and rather note that the one-body Hamiltonian 

Hi = -e^bl-bp + bl-h) 

m 

= ei Cm(e) , (58) 

gives the energy spectrum of a six-dimensional spherical oscillator 

Ei_(n) = em. (59) 

Here Ci[/(6) denotes the linear Casimir operator of U(6). The label n represents the number of oscillator 
quanta n = n p + ri\ = 0, 1, . . . , N. This special case is called the C/(6) limit. In Fig. || we show the 
structure of the spectrum of the spherical harmonic oscillator with U(6) symmetry. For three identical 
clusters, the wave functions must transform as the symmetric representation t = S of the permutation 
group £3. The levels are grouped into oscillator shells characterized by n. In Fig. ||, the levels belonging 
to an oscillator shell n are further classified by a — n, n — 2, . . . , 1 or for n odd or even. The quantum 
number a labels the representations of 5*0(6), a subgroup of f (6). The ground state has n = and 
L p = + . We note, that the n = a = 1 shell is absent, since it does not contain a symmetric state with 
t = S. The two-phonon multiplet n — 2 consists of the states L p = 2 + with a = 2 and + with a = 0. 
The degeneracy of the harmonic oscillator shells can be split by adding invariants of subgroups of U(6) 

For the six-dimensional spherical oscillator, the number of oscillator quanta n is a good quantum 
number. However, when vq in Eq. (p5[), the oscillator shells with An = ±2 are mixed, and the 
eigenfunctions are spread over many different oscillator shells. A dynamic symmetry that involves the 
mixing between oscillator shells, is provided by the reduction 

17(7) D SO(7) D SO(6) D ■ ■ ■ (60) 
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Again the dots indicate further subalgebras needed to classify the states uniquely. They can be taken as 
discussed previously. We consider now an algebraic Hamiltonian of the form 



Ho 



N(N + b)-(D p -D p + D x -D x ) 
N(N + 5) - C 2 so(7) + C 2 so(6) 

( S V - &+ . 6 t _ b \ . b [) (gg -h p -b p -h- l x ) + a 



2SO(6) 



where N is the number operator of Eq. (pq) and D denote the dipole operators 



D p.ra 
D\.m 



= (6+xs- 



fit x Mm 



(61) 



(62) 



The operators 0250(7) an d 2 so(6) represent the quadratic Casimir operators of SO (7) and 50(6), 
respectively. The energy spectrum in this case, called the £0(7) limit, is given by the eigenvalues of the 
Casimir operators (C 2 so(7)) = w ( w + 5) and (0 2 so(6)) = cr ( cr + 4) 



E 2 (uj,a) = k[{N -uj)(N + uj + 5) + cr((T + 4)] . 



(63) 



Here u> — N, N — 2, . . . , 1 or for N odd or even, respectively, labels the symmetric representations of 
50(7), and a — 0,2,3, ...,cj those of 50(6) (note that a = 1 is missing, since it doesnot contain a 
symmetric state) . The correspondence with the energy spectrum of the deformed oscillator of Eq. ( |33| ) 
can be seen by introducing a vibrational quantum number v — (N — u>)/2. The energy formula reduces 
to 



E 2 (v,a) = AnNvll- 



2v-5 
2N 



+ KCr(<7 + 4) 



(64) 



which to leading order in N, is linear in the vibrational quantum number. In Fig. |^ we show the spectrum 
of the deformed oscillator with 50(7) symmetry. The states are now ordered in bands labeled by ui, rather 
than in harmonic oscillator shells. Although the size of the model space, and hence the total number of 
states, is the same as for the harmonic oscillator, the ordering and classification of the states is different. 
For example, in the U(6) limit all states are vibrational, whereas the SO (7) limit gives rise to a rotational- 
vibrational spectrum, characterized by a series of rotational bands which are labeled by w, or equivalcntly 
by the vibrational quantum number v — (N — u>)/2 = 0, 1, . . . . 



4.2 Oblate symmetric top 

This situation does not correspond to a dynamic symmetry of the Hamiltonian H. However it is relatively 
simple to construct an algebraic Hamiltonian that produces the spectrum of an oblate symmetric top. 
This Hamiltonian was suggested in our preliminary account of this work H . It can be written as 



6 (sV 



b\)(~ss- 



bp ■ b p - b\ ■ b\) 



(bl ■ b^ p - b{ ■ b\)(b p ■b p -b x - b x ) + 4(6+ • b[)(bx ■ b p ) 



+6 

+2 Kl (bl xb p + b{x 6 A )« • (6+ x bp + 6+ x b x ) 



-3k 2 (6+ x b x 



6+ x S p )C) 



(b 



X x bp 



bl x b x ) 



(i) 

(0) 



(65) 



Since this Hamiltonian does not correspond to a dynamic symmetry, its spectrum must be obtained by 
diagonalizing it in an appropriate basis. A computer program has been written for this purpose by one of 
us (RB) and is available upon request ]l9| . An approximate energy formula can be obtained by making 
use of the method of intrinsic or coherent states as discussed in the following section. 
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5 Geometric analysis 



The geometric properties of the general algebraic Hamiltonian H of Eq. (|5^) can be studied by using 
intrinsic or coherent states |2(J ■ The ground state coherent state for the problem at hand can be written 
as 



\N;d p ,d\) 



1 



'JV! 



(bt) N \0) , 



with 



bt = \ II — dp ■ a * — a\ ■ a£ s ' + n 



p ■ bt +d x -b{ 



(66) 



(67) 



The condensate boson Eq. (|67]) is parametrized in terms of two vectors, a p and d\, which can be trans- 
formed to intrinsic coordinates q p , 9 P , 4> p , and q\, 9\, <p\, and their conjugate momenta pit] 



(68) 



with 



Pk,i 

Pk,0 

0k,-i 



[-Pfa/smOk - ipg k }/q k V2 
qk + ipk 

[-P4> k I sin 9 k +ipe k ]/qkV2 



(69) 



with k = p, A. The classical limit of a normal ordered operator can be defined as the coherent state 
expectation value of the operator divided by the number of bosons N. 

5.1 Dynamic symmetries 

We start by analyzing the two dynamic symmetries discussed in the previous section. For the U(6) limit, 
the classical limit of the Hamiltonian is given by 



Hi, 



1 

N 



(AT; dp, d\ | : Hi : | N; a p , a\) 



ei (dp ■ a* + d\ ■ a£) 

ei 5 (f P + $ + L 2 P /1 2 P + Pi + «A + L\/qi) , 



(70) 



where L 2 is the angular momentum in polar coordinates 



f 2 - 

Liu — 



(71) 



with k = p, A. Eq. (|7C| ) confirms the interpretation of the U(6) limit in terms of a six-dimensional 
spherical oscillator. 

For the 50(7) limit, we find 



H 



2,cl 



k(N- 1) 



if 



d 



L 2 L 2 

q 2 p 2 + (1 - q 2 ) 2 - ^ - 4cot2 X ^ + + — ^~ 

a\ ox sin x cos X 



Here we have made a change to hyperspherical variables 

q p = qsinx , q\ = qcosx 



(72) 



(73) 
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and their conjugate momenta. In the limit of small oscillations around the minimum of the potential 
energy surface q = 1 + Ag we find to leading order in N a one-dimensional harmonic oscillator in the 
hyperradius q 

H 2 , c i -> kN (p 2 + A(Aq) 2 ) , (74) 

with frequency AkN. To leading order in N, the frequency coincides with that of Eq. (|64|). This analysis 
shown the connection between the SO (7) dynamical symmetry and the displaced oscillator. 



5.2 Oblate symmetric top 



For the oblate symmetric top, the method of intrinsic states becomes particularly useful, since this case 
does not correpond to a dynamic symmetry. We first consider the vibrational part of the Hamiltonian of 
Eq. (|6i) 



H. 



3,vib 



= & (R 2 sV - b\ ■ b\ - b{ ■ b{) (R 2 ss -b p -b p -b x -b 



H2 



(K -bl~b\- b{) (b p ■b p -b x -b x)+4 (ftt . b \ ) (b x ■ b p ) 



(75) 



For R 2 = 0, this Hamiltonian has U(7) D U(6) symmetry and corresponds to a spherical vibrator, whereas 
for R 2 = 1 and £2 = it has U(7) D 50(7) symmetry and corresponds to a deformed oscillator. Here we 
analyze the general case with R 2 ^ and £1, £2 > 0. The classical limit of the Hamiltonian of Eq. ( [75] ) 
has a complicated structure. The equilibrium configuration of the potential energy surface is given by 



<7u 



\/2R 2 / (1 + R 2 ) , 



Xo 



tt/4 , 



Co = vr/4 , 



(76) 



where the relative angle between a p and d\ is denoted by 2£. In the limit of small oscillations around 
the minimum q = qa + Ag, x = Xo + A% and ( = Co + AC, the intrinsic degrees of freedom q, x an d C 
decouple and become harmonic. To leading order in N we find 



H. 



3,cl 



2R 2 



l + R 2 



p 2 + 2R 2 (l + R 2 )(AqY 
AR 4 



Pi 



AR 4 



r(AC)' 



(77) 



(l + i? 2 ) 2V ^' ^ (1+i? 2 ) 2 

Standard quantization of the harmonic oscillator yields the vibrational energy spectrum of an oblate top 



E: 



with frequencies 



3,vib 



4NR £1 



+ i) + UJ 2 (V2 + 1) , 



ANR 2 



lO-2 



(78) 



(79) 



in agreement with the results obtained in a normal mode analysis |15j . Here v\ and v 2 represent vibrational 
quantum numbers for a one- and two-dimensional harmonic oscillator, respectively. 
We consider next the rotational part of the Hamiltonian which can be rewritten as 



ff 3 ,rot = 2«i (6j xb p + b{ x& A )W-(6t xb p + b\ x& A )W 
+3^2 (bl xb x -b{x b p )W . (fit x l p _ b t x g A ) 



(0) 



(80) 



Both terms commute with the general S*3 invariant Hamiltonian of Eq. (J5^) and hence correspond to 
exact symmetries. The eigenvalues are given by 



E: 



3, rot 



= Hi L{L 



k 2 m 2 F 



n 2 (K T 2£ 2 ) 2 



(81) 
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Here we have used, that for the oblate top the quantum number mp is related to the projection K of the 
angular momentum on the symmetry-axis and £2 ]22] | . The last term contains the effects of the Coriolis 
force which gives rise to a 8k2-£^2 splitting of the rotational levels. 

The total (rotation-vibration) wave functions of a rigid triangular configuration can be written as 

\(vi,v£);K,L?,M). (82) 

For three identical particles, the wave functions for rigid configurations must transform as the symmetric 
representations A\ of T>^ ~ S3. This imposes some conditions on the allowed values of the angular 
momenta. For vibrational bands with £2 = and 1 , the allowed values of the angular momenta are 

(vi,-U2 2=0 ) : K = 3n n= 0,1,2,... 

L = 0,2,4,... iorK = 

L = K,K + l,K + 2,... for K ^0 (83) 
{vi,v{ 2=1 ) : K = 3n + l,3n + 2 n = 0,l,2,... 

L = K,K + l,K + 2,... 

The parity is P — (— ) K . The vibrational band (1, 0°) has the same angular momenta L p = + , 2 + , 3~~, 4 ± , . 
as the ground state band (0, 0°), while the angular momentum content of the doubly degenerate vibration 
(0, l 1 ) is given by L p = 1~, 2^, 3^, .... Combining all results one can write an approximate expression 
for the energy eigenvalues fo the oblate symmetric top 

E{v x ,vf,K,L,M) = £h+ Cl ; 1 (t, 1 + i)+ W2 (w 2 + i) 

+/ci L(L + 1) + k 2 (Kt 2£ 2 ) 2 ■ (84) 

In Fig. ^ we show the structure of the spectrum of the oblate top according to the approximate energy 
formula of Eq. (Q) . The energy spectrum consists of a series of rotational bands labeled by (v\ , v^ 2 ) and 
K. The degeneracy between states with different values of K can be split by the last term in Eq. (p4). 



6 Transition probabilities 

In order to calculate transition form factors in the algebraic cluster model one has to express the transition 
operator in terms of the algebraic operators. In the large N limit, the U(7) dipole operators D p and D\ 
of Eq. (|62] ) correspond to the Jacobi coordinates plf. The matrix elements of exp(— iqy / 2/3\ z ) can be 
obtained algebraically by making the replacement [ |15| 

^/2/3\ z -» PDx,z/X D , (85) 

where (3 represents the scale of the coordinate and Xjj is given by the reduced matrix element of the 
dipole operator |^5|. In summary, the transition form factors can be expressed in the ACM as 

T{i^f-q) = (jf,Lf,M I f I 7;, Li,M) , (86) 

with 

For the general ACM Hamiltonian of Eq. (^5|), the form factors cannot be obtained in closed analytic 
form, but have to be calculated numerically. Hereto, a computer program has been developed p9[ , in 
which the form factors are obtained exactly by using the symmetry properties of the transition operator 
of Eq. @ (see Appendix D of @). 
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6.1 Dynamic symmetries 



When a dynamic symmetry occurs the matrix elements of the transition operators can be obtained in 
explicit analytic form. The general procedure to derive the transition form factors in the U(6) and SO (7) 
limits has been discussed in |^3|, . 

In the U (6) limit, the elastic form factor is given by [ fl5"[ 

T(0+^0+;q) = {cos e) N ^ e- q2 132/6 . (88) 
In the last step we have taken the large N limit, such that eXjj — e%/3iV — — q(3 remains finite. For the 



extended charge distribution, the form factor of Eq. (88) is multiplied by a Gaussian exp(— q 2 /4a). The 
coefficients (3 and a cannot be determined independently, since both have the same dependence on q 2 . 
The charge radius can be used to fix one particular combination 



(rY /2 = \I^+P 2 - (89) 

In Table || we show the results for some transition form factors in the U(6) limit of the ACM. In the large 
N limit, they reduce to those obtained for the harmonic oscillator of Table |l|. All form factors exhibit 
an exponential fall-off with momentum transfer. In Fig. ^| we show the elastic form factor in the large N 
limit (solid line). The coefficients (3 and a were taken to give the charge radius of 12 C: (r 2 ) 1 / 2 = 2.468 
fm p4[ . The corresponding charge distribution shows an exponential decay with r (solid line in Fig. ^). 

For the 5*0(7) symmetry, the transition form factors can be expressed in terms of Gegenbauer poly- 
nomials or, equivalently, hypergeometric funcions. For the elastic form factor we have [ |T7j] 

H0t^0t;q) = ^Cffcose) 

N N + 5 2 4J 2 (g/3V2) 
- 2F 1 (- — ,—^— ,3; sin e) -» ^2 ■ ( 90 ) 



The normalization constant is given by Xr> = y/N(N + 5)/2. Again, in the last step we have taken the 
large N limit, such that eXjj = —qf3 remains finite. The SO(7) transition form factors are presented in 
Table |[ In the large N limit, they reduce to those obtained for the deformed oscillator of Table |[ All 
form factors exhibit an oscillatory behavior. The transition form factor to the first excited + state which 
belongs to a vibrational excitation v = (N — oS)/2 = 1 vanishes identically, since the dipole operator is a 
generator of SO(7). In Fig. | we show the elastic form factor for N — > oo (dashed line). The values of (3 
and a are taken to reproduce the minimum in the elastic form factor of 12 C at q 2 = 3.4 fm~ 2 , and the 
12 C charge radius of 2.468 fm. The corresponding charge distribution for the deformed oscillator is given 
by the dashed line in Fig. ^. 



6.2 Oblate symmetric top 

For the oblate symmetric top, the form factors can only be obtained in closed form in the large N limit 
pl| (see Table [l]). As an example, the elastic form factor has been obtained as [|l5], |l8| 

^(0+^0+;<z) jo{q0). (91) 

In this case, the normalization constant is given by X D = NRy/2/(l + R 2 ). As can be seen from the 
last column of Table [j], all form factors show an oscillatory behavior. In Fig. ^ we show the elastic form 
factor in the large N limit (dotted line) . The coefficients (3 and a have been determined to give the same 
value of the minimum and the charge radius as for the deformed oscillator. The charge distribution for 
the oblate top is given by Eq. ([50]) , and is represented by the dotted line in Fig. ||. 
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The transition form factors for vibrational excitations are proportional to the corresponding intrinsic 
transition matrix elements xi an( i Xi f° r each type of vibration v\ and V2, respectively H 



l-R 2 



Xi 



2RVN 



X2 = 7=- • (92) 

Ry/2N 

In Fig. ^ we show the inelastic form factor 1^(0^ — > 0^";g)| 2 for different values of R 2 . As before, the 
coefficients (3 and a are determined by the first minimum of the elastic form factor and the charge radius. 
Whereas the elastic form factor does not depend on R 2 , the inelastic form factor is very sensitive to the 
value of i? 2 , especially with respect to the position of its minimum. 



7 The nucleus 12 C 

The structure of 12 C has been extensively investigated since the early days of Nuclear Physics. Compre- 
hensive calculations of 12 C have been done within the framework of the shell model, starting with the 
early calculation of Cohen and Kurath p5| in the p-shell and ending with the more recent calculations 
with large shell mixing. 12 C has also been the battleground for investigations within the framework of the 
microscopic cluster model. Here we investigate the extent to which properties of the low-lying spectrum 
of 12 C can be described in terms of the macroscopic cluster model. The reason why a detailed analysis 
is now feasible is two-fold. First, we have a set of explicit analytic formulas for energies, electromagnetic 
transition rates and form factors for the oblate top that can be easily be compared with experiments. 
Second and foremost is that the development of the algebraic cluster model (ACM) allows us to study in 
a straighforward way complex situations, such as those that are not rigid, but correspond for example to 
situations which are intermediate between two limiting cases. 

The differences between the three special solutions of the ACM, the U(6) limit, the SO (7) limit and 
the oblate top, are most pronounced in the transition form factors. The elastic form factor of 12 C shows 
a minimum at q 2 = 3.4 fm~ 2 and no further structure up to about q 2 ~ 14 fm~ 2 . Also the inelastic 
form factor for the transition — > 0^ shows a clear minimum at ~ q 2 — 4.5 fm~ 2 . The analysis of the 
form factors that was presented in the previous section shows clearly that only the oblate top scenario 
can account for the features of the empirical form factors of 12 C. In the U(6) limit the form factors 
drop off exponentially and have no minima, whereas in the SO(7) limit the inelastic form factor vanishes 
identically. 

Therefore, in this section we analyze the spectroscopy of 12 C in the oblate top limit of the ACM. 



7.1 Energies 



The Hamiltonian for the oblate top is given by Eq. (65). The coefficients £1, £2, K i and K2 are determined 
in a fit to the excitation energies of 12 C (see Table|). The number of bosons is taken to be N = 10. 
In Fig. H we show a comparison between the experimental data and the calculated states of the oblate 
top with energies < 15 MeV. One can clearly identify in the experimental spectrum the states + , 2 + , 
3 _ , 4 + of the ground rotational band, the state + of the stretching vibration and the state 1~ of the 
bending vibration. However, the 3 _ state does not fall at the location expected from a rigid extended 
configuration with moments of inertia given by Eq. (0), but at a higher excitation energy. Since the 
3 _ state has K = 3 the deviation from the simple rotational formula indicates large rotation-vibration 
interactions. This conclusion is strengthened by the low location of the stretching vibration. In contrast 
with molecules, the vibrational and rotational frequencies in nuclei appear to be comparable. If one takes 
as a measure of the ratio E V ib/E Iot the quantity E Q +/E 2 +, this quantity is 1.7 in 12 C and ~ 10 in a 
typical molecule with D^h symmetry (for example, ozone O3). 
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7.2 Form factors and transition rates 



Form factors for electron scattering on 12 C have been measured long ago. For the extended charge 
distribution of Eq. the form factors are obtained by multiplying the oblate top form factors by an 
exponential factor exp(— q 2 /4a) (compare Eqs. ( po|) and ( |24|)) . The coefficient f3 is determined from the 
first minimum in the elastic form factor at q 2 = 3.4 fm~ 2 ]24| to be j3 = 1.74 fm, and subsequently the 
coefficient a is determined from the charge radius of 12 C, (r 2 ) 1 / 2 = 2.468± 0.012 fm @, to be a = 0.52 
fm~ 2 . If the alpha-particles would not be effected by the presence of the others, the coefficient a should 
be the same as for free particles. The value of a is slightly different from the one used in ||, where it 
was determined from the maximum in the elastic form factor. 

In Fig. H we show a comparison between experimental and theoretical form factors calculated exactly 
using the Hamiltonian of Eq. ( |65| ) . The value of the parameters of the Hamiltonian and the transition 
operator is given in Table |]. We note that the transition form factors are calculated with the same values 
of coefficients a and /3 as determined from the elastic form factor. Whereas for the rotational excitations 
the calculations seem to give a good description of the experimental data, the situation is different for 
vibrational excitations. Although the q dependence of the vibrational form factors is consistent with 
experiment, indicating that the 0^" and lj~ excited states could indeed be the vibrational excitations of a 
three-alpha configuration with T>3h symmetry, their observed large strength implies a very large mixing 
with other configurations. The strengths of these excitations are proportional to the corresponding 
intrinsic transition matrix elements \i an< i Xi f° r each type of vibration v\ and respectively ||. 
According to Eq. (p2|), the coefficients xi an d \2 g° to zero in the large N limit (rigid configurations). 
For finite values of N, the vibrational transition form factors are different from zero. The solid lines in 
Figs. H were obtained in calculations for a finite value of N = 10. 

From the form factors one can extract other electromagnetic properties of interest. The ground state 
charge distribution can be obtained by taking the Fourier transform of the elastic form factor according 
to Eq. (|49|). In Fig. [l(] we compare the charge distribution of 12 C obtained from the calculated elastic 
form factor of Fig. ^ and extracted in a Fourier-Bessel analysis of the experimental data 0|. 

The B(EL) values correspond to the long wavelength limit of the form factors (see Eq. (5^)). The 
values extracted from the fit to the form factors are shown in Table 0. For the monopole strength, 
instead of B(E0) it is customary to give the matrix element M(E0). While the B{E2\2\ — » Q+) and 
B(E3;3i -> 0^) values are in reasonable agreement, both the B(E2;0% -> 2^) and M{EQ\G$ -> Of) 
deviate by an order of magnitude. This again may indicate that the nature of the 0^ is not that of a 
vibrational mode. 



8 Conclusions 

In this article, we have presented a detailed analysis of a configuration composed of three alpha particles 
at the vertices of an equilateral triangle. We have done so in a classical, quantum and quantum algebraic 
description. The formulas that have been derived have been used to study the ground state configura- 
tion of 12 C. In addition, we have introduced a new method to study clustering based on the algebraic 
quantization of the Jacobi variables. The latter method is very general and it allows to describe in a 
relatively simple fashion not only three particles at the vertices of a triangle (oblate top) but also other 
situations such as six-dimensional vibrations and rotational- vibrational spectra. 

Our analysis of the experimental data for 12 C appears to indicate that the triangular configuration 
describes the data (energy levels and form factors) reasonably well for the rotational band Of, 2f, 3^, 
4^, although with large rotation- vibration interactions. The situation is different for the vibrational 
excitations 0^ and 1]~. Here the shape of the form factors is well reproduced but its magnitude is not. 

We bring attention to the quantum algebraic description in terms of the algebra of U(7) that is capable 
of describing not only rigid like structures but also floppy structures. We point out that this method is 
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very general and can describe other situations, as, for example, three alpha particles on a line. These 
configurations cannot describe the low- lying experimental spectrum of 12 C since for three alpha particles 
on a line there are (in the ground rotational band) no negative parity states. They may occur as excited 
configurations. Their properties will be presented in a forthcoming publication. The crucial problem 
in clustering in nuclei is to understand what type of configurations are present, if any, and to provide 
unambiguous experimental evidence for these configurations. The algebraic method gives a way in which 
all calculations can be performed easily. 

Finally, the general algebraic framework based on the spectrum generating algebra U (7) also allows 
to describe clustering phenomena for three different particles J34|. This description is relevant to giant 
trinuclear molecules in ternary cold fission |35], |36| . 

Acknowledgements 

It is a pleasure to thank Aurora Tumino for her help in calculating the form factors of 12 C. This work 
was supported in part by CONACyT under project 32416-E, and by DPAGA under project IN106400, 
and by D.O.E. Grant DE-FG02-91ER40608. 

References 

[1] JA. Wheeler, Phys. Rev. 52 (1937), 1083. 

[2] K. Wildermuth and Th. Kanellopoulos, Nucl. Phys. 7 (1958), 150. 

[3] K. Ikeda, H. Horiuchi and S. Sato, Suppl. Progr. Theor. Phys. 68 (1980), 1. 

[4] Y. Fujiwara, H. Horiuchi, K. Ikeda, M. Kamimura, K. Kato, Y. Suzuki and E. Uegaki, Suppl. Progr. 
Theor. Phys. 68 (1980), 29. 

[5] H. Horiuchi, Prog. Theor. Phys. 51 (1974), 1266; Prog. Theor. Phys. 53 (1975), 447. 

[6] F. Iachello, in "Lie Algebras, Cohomologies and New Applications of Quantum Mechanics" (N. 
Kamran and P. Olver, Eds.), Contemporary Mathematics, AMS, Vol. 160, p. 151, 1994. 

[7] F. Iachello, Phys. Rev. C 23 (1981), 2778. 

[8] R. Bijker and F. Iachello, Phys. Rev. C61 (2000), 067305. 

[9] A. Tumino, R. Bijker and F. Iachello, in "Proceedings of the 7th International Conference on Clus- 
tering Aspects of Nuclear Structure and Dynamics" (M. Korolija, Z. Basrak and R. Caplar, Eds.), 
p. 271, World Scientific, Singapore, 2000. 

[10] G. Herzberg, "Molecular Spectra and Molecular Structure II: Infrared and Raman Spectra of Poly- 
atomic Molecules", Van Nostrand, New York, 1950. 

[11] J.D. Jackson, "Classical Electrodynamics", Wiley, New York, 1962. 

[12] F.T. Smith, Phys. Rev. 120 (1960), 1058. 

[13] J.L. Ballot and M. Fabre de la Ripelle, Ann. Phys. (N.Y.) 127 (1980), 62. 

[14] G. De Franceschi, F. Palumbo and N. Lo Iudicc, Phys. Rev. C 29 (1984), 1496. 

[15] R. Bijker, F. Iachello and A. Leviatan, Ann. Phys. (N.Y.) 236 (1994), 69. 



18 



[16] L.A. Copley, G. Karl and E. Obryk, Nucl. Phys. B 13 (1969), 303; 
F.E. Close and Z. Li, Phys. Rev. D 42 (1990), 2194. 

[17] R. Bijkcr and A. Leviatan, Rev. Mex. Fi's. 44 S2 (1998), 15. 

[18] E.V. Inopin and B.I. Tishchenko, Soviet Physics, JEPT 11 (1960), 840. 

[19] R. Bijker, Computer program ACM, unpublished. 

[20] J.N. Ginocchio and M. Kirson, Phys. Rev. Lett. 44 (1980), 1744; 

A.E.L. Dieperink, O. Scholten and F. Iachcllo, Phys. Rev. Lett. 44 (1980), 1747. 

[21] O.S. van Roosmalen and A.E.L. Dieperink, Ann. Phys. (N.Y.) 139 (1982), 198. 

[22] R. Bijker, A.E.L. Dieperink and A. Leviatan, Phys. Rev. A 52 (1995), 2786. 

[23] R. Bijker and J.N. Ginocchio, Phys. Rev. C 45 (1992), 3030. 

[24] W. Reuter, G. Fricke, K. Merle and H. Miska, Phys. Rev. C 26 (1982), 806. 

[25] S. Cohen and D. Kurath, Nucl. Phys. 73 (1965), 1. 

[26] F. Ajzenberg-Selove, Nucl. Phys. A 506 (1990), 1. 

[27] I. Sick and J.S. McCarthy, Nucl. Phys. A bf 150 (1970), 631. 

[28] H.L. Crannell and T.A. Griffy, Phys. Rev. 136 (1964), B1580. 

[29] H. Crannell, Phys. Rev. 148 (1966), 1107. 

[30] H. Crannell, J.T. O'Brien and D.I. Stober, "Int. Conf. on Nuclear Physics with Electromagnetic 
Interactions" , Mainz, 1979. 

[31] P. Strehl and Th.H. Schucan, Phys. Lett. 27B (1968), 641. 

[32] Y. Torizuka, M. Oyamada, K. Nakahara, K. Sugiyama, Y. Kojima, T. Terasawa, K. Itoh, A. Yam- 
aguchi and M. Kimura, Phys. Rev. Lett. 22 (1969), 544. 

[33] A. Nakada, Y. Torizuka and Y. Horikawa, Phys. Rev. Lett. 27 (1971), 745 and 1102. 

[34] R. Bijker and A. Leviatan, Few-Body Systems 25 (1998), 89. 

[35] §. Misicu, P.O. Hess and W. Greiner, Phys. Rev. C 63 (2001), 054308. 

[36] R. Bijker, P.O. Hess and §. Mi§icu, Heavy Ion Physics 13 (2001), 89; 
P.O. Hess, R. Bijker and §. Mi§icu, Rev. Mex. Fis. 47 S2 (2001), 52. 



19 



Table 1: Transition form factors .F(0^ — > Lf ; q) for the harmonic oscillator (ho), the deformed oscillator 
(do) and the oblate top (ot). For each case, a scale parameter /3 was introduced to give the same result 
for the charge radius. 



Lf 


ho 


JT(0+ ^ if ; <?) 
do 


ot 


ot 


e -9 2 /3 2 /6 


4J 2 (qpV2)/q 2 P 2 


Jo (9/3) 


ot 









ir 




z8V3J 5 ((z/3V2)/ g 2 /3 2 


-iX2lV3q(3j 2 (qf3) 




^0 2 «2 -g 2 /3 2 /6 


-8V2J 4 (g/3V2)/(7 2 /3 2 






i 3^3 - 9 2 /3 2 /6 


z8V2J 5 ((z/3V2)/ g 2 /3 2 




4f 




^Lj 6 (q(3V2)/q 2 p 2 





Table 2: Transition form factors ^"(0^ — > if ; g) for the C7 (6) limit of the algebraic cluster model. The 
coefficient e is given by e = —q(3/y/3N. 



Lf HOt-Lf-q) 



Ot (cose) w 

2 + (cos e)^ sine) 2 

l i \f 40(N-3y. (cos e)^' 3 (» sine) 3 

2 f yi^ (cose)JV " 2(isine)2 

3 r \/ 60(^-3)! ( C0S £ ) Ar_3 (* Sin £ ) 3 

4 t \/ 280(^-4)! ( COS £ ) N ~ 4 (* Sin £ ) 4 
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Table 3: Transition form factors ^(Of — > if ; g) for the SO(7) limit of the algebraic cluster model. The 
coefficient e is given by e = — y / N(N + 5)/2. 



^(0+^Lf;g) 



u i (iv+4)! °jv (cose) 



0+ 



,- / l89jV!(jV-3)!10r r . ■ x 3 ^(11/2), x 



„+ / 14QjV!(jV-2)!8! ■ x 2 r (9/2) / x 



o- / 126jV!(jV-3)!ToT / ■ ■ x 3 r (ll/2) / x 

d l Y (Ar+4)!(W+7)! x* &me J °Ar-3 l cos e J 



4+ / 297jV!(jV-4)!T2T / ■ ■ x 4 r (13/2) / x 

4 1 V (iV+4)!(JV+8)! l z&me i °W-4 (CO^eJ 



Table 4: CoefBcients used in the calculation of the energy spectrum and the transition form factors of 
12 C. The number of bosons is N = 10. 



a 


0.1721 


MeV 


6 


0.2745 


MeV 


«i 


0.7068 


MeV 


K2 


0.1276 


MeV 


i? 2 


1.40 




/? 


1.74 


fm 


a 


0.52 


fm- 2 
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Table 5: Comparison between calculated and measured B(EL) values. 



Th. 



Exp. 



Rcf. 



B(E2;2+ -> 


of) 


8.4 


7.6 ±0.4 


e 2 fm 4 


5(£3;3J~ -> 


of) 


44 


103 ±17 


e 2 fm 6 


B(E4;4f -> 


of) 


73 




e 2 fm 8 


B(£2;0f -> 


2f) 


1.3 


13.1 ±1.8 


e 2 fm 4 


M(£0;0f - 

( r 2)l/2 


•of) 


0.4 


5.5 ±0.2 


fm 2 




2.468 


2.468 ±0.12 


fm 



m 

[31] 

[24] 




2 ~ ~ 3 

Figure 1: Geometry of a three-body system. 
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Figure 2: Schematic spectrum of the harmonic oscillator with U(6) D SO (6) symmetry. The number of 
bosons is TV = 4. All states are symmetric under 53. 
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Figure 3: Schematic spectrum of a deformed oscillator with SO (7) symmetry. The number of bosons is 
N = 4. All states are symmetric under S3. 
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Figure 4: Schematic spectrum of an oblate symmetric top. All states are symmetric under S3. 
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Figure 5: Comparison of the elastic form factor |.F(0^~ — > Of; q)\ 2 , calculated for N — > oo (a) in U{6) with 
= 1.70 fm and a = 0.47 fin" 2 (solid line), (b) in 50(7) with = 1.97 fm and a = 0.67 fm" 2 (dashed 
line), (c) in the oblate top with = 1.70 fm and a = 0.47 fm -2 (dotted line). 
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Figure 6: Ground state charge distribution calculated for N — > oo (a) in U(6) with f3 = 1.70 fm and 
a = 0.47 far 2 (solid line), (b) in 5*0(7) with (3 = 1.97 fm and a = 0.67 far 2 (dashed line), (c) in the 
oblate top with (3 = 1.70 fa and a = 0.47 far 2 (dotted line). 
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Figure 8: Comparison between the low- lying experimental spectrum of 12 C |26j and that calculated with 
the oblate top Hamiltonian of Eq. ( |65| ) with N — 10. The parameter values are given in Table [|. States 
with uncertain spin-parity assignment are in parentheses. 
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Figure 9: Comparison between the experimental form factors |.F(0* — > Lf ; q)\ 2 of 12 C for the final states 
(a) Lf = Of (elastic), (b) Lf = 2+, (c) Lf = 3jf , (d) Lf = 4+, (e) Lf = 0+, and (f) Lf = 1^ and those 
obtained for the oblate top with N = 10. The parameter values are given in Table |I| The experimental 
data are taken from il and K^-il. 
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Figure 10: Comparison between the ground state charge distribution of 12 C extracted from the experi- 
mental elastic form factor (dashed lines) [p4| and that calculated for the oblate top with N = 10 (solid 
line) . The parameter values are given in Table ||. 
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